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Abstract 



We present a lattice Monte Carlo simulation for the evaluation of the spec- 

K^ i trum of a confining Yang-Mills theory without Goldstone boson. We show that 

this theory is a very good candidate for describing composite weak bosons. In 

^ ■ order to perform the spectrum analysis we have used standard lattice QCD 

Monte Carlo methods. We have also developed an efficient method to evaluate 

the mass of the pseudoscalar isosinglet which is present in our theory. 



1 Introduction 

The Standard Model |[l|] (SM) describes the strong, weak and electromagnetic inter- 
actions by a gauge theory based on the group Gsm = SU{?>) x SU(2) x U{1) which is 
broken by the Higgs mechanism to SU{?)) x f/(l). The theory is essentially determined 
once the matter fields and their transformation under the local gauge transformations 
of Gsm are specified. The matter fields (leptons and quarks) and the Higgs boson 
are considered to be elementary. They interact with each other by the exchange of 
gauge bosons which are also considered to be elementary. The structure of the SM 
has been phenomenologically confirmed to high accuracy. 

In spite of the beautiful corroboration of the SM by experiments a natural ques- 
tions arises: How elementary are the leptons, the quarks, the Higgs bosons and the 
gauge bosons? The idea that the SM itself is an effective theory of another, more 
fundamental, where quarks, leptons and bosons are composites of more fundamental 
fields is almost as old as the SM itself. The idea of quark and lepton compositeness 
is motivated by the observed connection between quarks and leptons, by the gen- 
eration puzzle and by the existence of too many parameters in the SM. The Higgs 
compositeness is motivated by the fine tuning problem. The W and Z compositeness 
is motivated by their relation to a composite Higgs and by the observation that all 
short-range interactions are residual interactions of a more fundamental long-range 
interaction. 

The constituents, the new fundamental fields, are supposed to carry a new internal 
quantum number (which we denote as hypercolor) and the quarks, leptons and bosons 
are hyper color less composite systems of them. The binding of the constituents due 
to hypercolor is viewed as an analogy to the color confinement mechanism of QCD. 
However, since the SM spectrum is different from the hadron spectrum, the hyper- 
color interaction has to be described by a strongly coupled Yang-Mills theory different 
from QCD. 

Several models treat the quarks, leptons and bosons as composite systems. Today a 
conspicuous number of theorems exist which have ruled out most of the existing mod- 
els and radically restricted the possibilities to construct realistic composite models. 
In principle there are two categories of models: 

• Three-fermion models: In these models the leptons and quarks have three 
fermion spin 1/2 constituents. Most of these models are phenomenologically 
ruled out by the Weingarten 0], Nussinov |^ and Witten constraints, 
by the the 't Hooft's anomaly-matching conditions 0, by the Weinberg- Witten 
theorem [0] or by the Vafa- Witten theorem |^ . Because most of the constraints 
dictated by these theorems can be avoided in a supersymmetric scenario of 
compositeness, some authors have proposed supersymmetric versions of these 
models & ITOi. 



Fermion-scalar models: In these models the leptons and quarks are made 
out of two constituents, a fermion spin 1/2 and a scalar. The W, Z and the 
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Higgs bosons can have also two constituents. Two particular examples are the 
"Strongly coupled SM" (SCSM) |lTl and the Yang-Mills theory without Gold- 
stone bosons [T^Q- 

The SCSM starts from the same lagrangian of the SM, however, with an un- 
broken gauge group SU{2) ® U{1). The aim of this theory is not to describe 
a model more fundamental than the SM, but to propose an alternative to the 
usual Higgs mechanism. The interesting feature arises that the W and Z are 
only the lowest lying states of a whole spectrum of vector-mesons. 
The Yang-Mills theory without Goldstone bosons is a strongly coupled non- 
abelian gauge theory which can have a low lying spectrum equal to the spec- 
trum of the SM. This is the model we want to study in this work. In particular 
we will concentrate on the weak gauge bosons sector comparing the low lying 
spectrum of this theory with the SM boson spectrum. 

To be precise this last model considers the photon to remain elementary and switched 
off. The weak gauge bosons W^ and Z^ then form a mass degenerate triplet. This 
model is a usual confining Yang-Mills theory with SU{2) local hypercolor gauge group, 
SU{2) global isospin group and generalized Majorana fermions in the fundamental 
representation of the local and global symmetry groups. We note that the general- 
ization of the Majorana fermions with non-trivial quantum numbers is possible only 
if the symmetry groups are real. 

To be viable a composite model of the weak bosons has to reproduce the known weak 
boson spectrum: the lightest bound states have to be the W-bosons and heavier bound 
states have to lie in an experimentally unexplored energy range. The only possibility 
to have a Yang-Mills theory which reproduces the weak boson spectrum is to choose 
the degrees of freedom in a way that they naturally avoid bound states lighter than 
the vector isotriplet of the theory which characterizes the W-boson triplet. This is 
possible if the unwanted light bound states which naturally show up as Goldstone 
bosons or pseudo Goldstone bosons in many models (like, for example, a pseudoscalar 
isomultiplet, which would be the pion analogue of QCD) are avoided. The choice of 
Majorana fermions in this model avoids the SUaC^) global chiral symmetry of the 
Yang- Mills Lagrangian because left- and right-handed degrees of freedom are not in- 
dependent. The axial current (which would generate the SUaC^) chiral symmetry) 
does not exist and it is not possible to have a breaking of SUAi'i) with the related 
low lying Goldstone bosons. In fact, the pseudoscalar isotriplet vanishes by the Pauli 
principle (it is a symmetric combination of Grassmann variables). 

Because of the strong coupling character of this theory, we need non-perturbative 
methods to make predictions. It is important that the fermion theory under discus- 
sion can be defined by a gauge invariant lattice regularization. A lattice regularization 



"'^In |I2] an unfortunate sign error has entered the calculations due to a wrong factor i in the 
fourth and fifth terms of eq. (3). To interpret the results correctly one has to interchange isosinglet 
<-> isovector in n%. I thank H.Schleret for a discussion of this point. 



a la Wilson [JT^I is possible because the choice of the isospin group SU{2) allows us to 
replace the Dirac mass term and the Dirac-type Wilson term by a hypercolor gauge 
invariant Majorana type expression. 

A strong coupling expansion analysis |jl^ of the spectrum of this theory has shown 
that the spin one isotriplet bound state (the right quantum number to represent the 
W-boson of the SM) could be the lightest state if the pseudoscalar isosinglet acquires 
a mass by the chiral anomaly in analogy to the r]' in QCD. 

In this work we calculate the spectrum of the lightest bound states by a quenched 
Monte Carlo simulation and we show that the vector isotriplet bound state of this 
theory is the lightest one. We have performed two different type of Monte Carlo 
simulations. In the first one we have calculated the masses of bound states without 
the contribution of the chiral anomaly. This simulation confirmed the results of the 
strong coupling expansion. In the second one we have developed an efficient method 
to evaluate the chiral anomaly contribution to the mass of the pseudoscalar isosinglet 
bound state. Its mass turned out to be heavier than the vector isotriplet mass. 

Our work is organized as follows: in section 2 we introduce the model in ques- 
tion and because this paper is not addressed only to lattice specialists in section 3 
we shortly present the technical methods of lattice spectroscopy that we will use. 
In section 4 we explain the method that we have developed to estimate the chiral 
anomaly contribution to the pseudoscalar isosinglet mass and in section 5 we analyse 
the results of the simulations. 

2 Confining Gauge Theories without Goldstone Bosons 

2.1 The Wilson action 

We consider a gauge theory whose fermion content is represented by a Weyl spinor 
F^^{x). Here a denotes the (undotted) spinor index (a = 1,2), A denotes the 
fundamental representation index of a global SU(2) isospin group {A = 1, 2) and a 
denotes the fundamental representation index of the local SU(2) hypercolor gauge 
group (a = 1,2). We introduce the generalized Majorana spinor ip starting from the 
Weyl spinors F and its conjugate F'^ 





and its adjoint 

V^(x) = (F^(x)Q,Ft(x)) = ¥.^(x)f "i ''] (2) 




The matrix Q represents the antisymmetric matrix in spin, hypercolor and isospin 
space, which correspond to the Kronecker product of ia2,iT2,iT2 (the antisymmetric 
matrices in spin, hypercolor and isospin space, respectively). Of course the fields ip 



and ip are not independent fields. The choice of the global isospin group SU{2) and 
of the local hypercolor group SU{2) allows us to write gauge invariant mass terms 
for the Majorana fermion fields 

V^V' = FQF + F^QF^ (3) 

Note that this choice is unique if one deals with Majorana fermions. Because of 
the existence of the mass term we can define the Yang-Mills action on the lattice in 
Euclidean space in the form of a Wilson action. 

S = /3j2TrUidp)-k Y^ ^ix)Tib)Uib)^Piy)-W,Pix)^P{x)^ 

pCA b={xy) X 

^ SGauge{U)-i'^^i^ (4) 



For further details we refer to ref. 14 . 



We will make repeated use of the fact that the fermion fields can be integrated out 
from the functional integral. Here, one has to be careful since ip and ip are not 
independent fields and the fermionic functional integral is represented in terms of (p. 
As a result we obtain. 

Z = J[d^][dU]exp{-S{ip,U)} = j[dU]^d^tM(jr)ew{-SGauge{U)] (5) 

2.2 CP eigenstates 

In weak interactions CP is a good quantum number but not C and P separately. 
Therefore, we perform a classification of the composite operators according to the 
CP eigenvalues. The CP transformation of the Weyl spinor F is defined by 

F^''{x)=iT2ia2F\xp) 

(where xp = {t, —x)) which fixes the phase. In the compact notation we can rewrite 
the CP transformation in the following form: 

^^^(a;) = T27o^(xp) 

V'^^(x) = ^(xp)7oT2 (6) 

We can build CP odd and CP even eigenstates^. In a confining Yang-Mills theory the 
physical states are hypercolor singlets. Because in our model the hypercolor group is 
SU{2) all physical states are described by bilinear forms in ip and ip. 

^Lorentz scalars S, vectors V^" and tensors Tt"" are CP even if CPS = S, CP V = V,, and 
CPTt^" = T^^ and are CP odd iiCPS ^ -S, CPVt' = -V^ and CPTi^" = -T^^, respectively. 



• Scalar CP eigenstates: The scalar CP even combination of scalars is 

S+{x) = ■ijj{x)ip{x) (7) 

The pseudoscalar CP odd combination is 

S^{x) = itlj{x)'y^ip{x) (8) 

• Vector CP eigenstates: The axial vector isosinglet CP even combination is 

\/^(x) = lv^(x)7^75^(x) (9) 

The isotriplet vector states is 

V^'^x) = ^ij{x)-f^'T^ij{x) (10) 

This state is CP even for 1=1,3 and CP odd for 1=2 and characterizes the 
W-boson triplet. 

• Tensor CP eigenstates: The tensors are 

B"''^ =ilj{x)a'"'T^iP{x) (11) 

and their dual 

5f ^ = ef"'''^'Bi^ (12) 

where a'^'^ = ^['y^,'^^]- B'^'^^ has CP properties opposite to the vector isotriplet 
V^ ones. Its dual B^"^ has CP properties opposite to B^^"^ . 

2.3 Correlators 

Masses are computed in lattice simulation from the asymptotic behavior in the 
Euclidean-time direction of static correlation functions. Static here means that we 
have projected the correlation functions at zero momentum transfer. In our case we 
need to consider only two point functions. A typical static bound state propagator 
can be written as 

Crit)=J2{{i,it,x)Ti,it,x)ymTi,iO)) (13) 

X 

where the F represent a matrix with hypercolor, isospin and spinor indices, for exam- 
ple r = 'j^T^ for the vector isotriplet. Here the contraction of the hypercolor, isospin 
and spinor indices is implicit. Contracting creation and annihilation operators into 
fermion propagators 

{T^Pix)^Piy)) = Gix,y) (14) 



we obtain for the correlation function 

Crit) = 5:((v^(x)rv^(x))tv;(i/)rv^(i/)) = (15) 

X 

= J2 {^^ {^G{x, 0)rG(0, x) - TG{x, 0)AT^AG{0, x)) - 

X 

- (Tr [TG{0, 0)] Tr [TG{x, x)] - Tr [aT^AG{0, 0)] Tr [aT^AG{x, x)] ) } 

where T^ indicates the transposition of F and A is the matrix A = Q(BQ- Siniphfying 
this expression we obtain the full propagator of the CP eigenstates. 

Cr{t) = 2j2{Tr [rG{x, 0)rG(0, x)] - Tr [rG(0, 0)] Tr [rG{x, x)]} (16) 

X 

To prove the formula (15) one has to contract properly all fermion fields in (13) 
rewriting them in terms of if fields and remembering that they are not Dirac fields. 
The disconnected part of eq. (16) corresponds to the chiral anomaly contribution to 
the propagators of the CP eigenstates. This contribution vanishes for all isomultiplet 
bound states. It yields an important contributions for the pseudoscalar isosinglet 
bound states S-. 

The anomaly contribution may be switched off by restricting to the sector of configu- 
rations of topological charge zero. We will make use of the fact that when the chiral 
anomaly contribution is switched off the pseudoscalar isosinglet bound state behaves 
like a massless Goldstone bosons and has a propagator of the form 

C^3(t) =2 5: Tr [75^(^,0)75^(0, x)] (17) 

X 

On the other hand, when we take into account the chiral anomaly contribution it is 
no more a massless Goldstone boson but it acquires a mass, like the rj' in QCD. Its 
propagator is of the form 

G,7(t) = 2 X: {Tr [75G(x, 0)75^(0, x)] - iNjTr [75^(0, 0)] Tr [^,G{x, x)] } (18) 

X 

Here Nj = 2 characterizes the SU{Nf) isospin group. Notice that we have introduced 
the parameters ^ in the propagator (18). This parameter is ^ = 1 if one computes 
with dynamical fermions. Unfortunately, present computational capabilities make 
necessary the quenched approximation and the parameter ^ can acquire a value dif- 
ferent from unity |]I3| . 



We will evaluate the mass once without and once with anomaly term. In the case 
without the chiral anomaly contribution the mass of the operator |^ S- will be used 
to identify the chiral limit. The propagator C^^{t) can be evaluated using stan- 
dard Monte Carlo techniques. For the case where the chiral anomaly contribution is 

^To distinguish the two cases we denote by S";^ the pseudoscalar isosinglet when the chiral 
anomaly is switched on. 



switched on we have used an improved version of the method of ref. [|r3| to evaluate 
the propagator C^"(t). 

Because we will evaluate static propagators we consider the operators being given 
by the following set of F's.Q 
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3 Lattice spectroscopy 

The general procedure in lattice spectroscopy is to first calculate the fermion prop- 
agator G{x, y) in the presence of an external gauge field and then to evaluate the 
hypercolor singlet propagators (16). The action (4) is quadratic in the fermion fields, 
these can be formally integrated out and the fermion propagator needed to evaluate 
the correlation functions (16) takes the form: 



Gix,y) = {T^ix)^Piy)) = - [d^][dU]Ti;ix)meM-Si^,U)} 



[dU] M{U)~l^detM{U) eM-SaaugeiU)} 



(19) 



Just as a simple integral can be evaluated as a limit of a sum, we can evaluate the path 
integral (19) by discretizing it on a lattice in an Euclidian four dimensional space. 
We put our model on a iV^ x A^^ lattice with spacing a. A gauge field configuration 
is defined as a set of A^ x A^ complex matrices^ defined on each oriented lattice bond 
h with some boundary condition. Fermion fields are defined as A^ x Nf x 4 complex 
vectors on the lattice points with some boundary condition. The lattice spacing a 
acts as ultra-violet cut-off and provides a regularization scheme necessary for any 
quantum field theory. The Wilson lattice action (4) provides a regularization which 

^Because of the relations d^V^ — C x S- (when the anomaly is switched off) and d^B^^"^ = 
C' X V'^^ (where C and C" are constants) the bound states F° and B^^^ are not independent from 5_ 
and, respectively, V^ . Therefore, after summing over all x in eq. (13) we expect that the correlators 
of V^ and B^^^ receive contributions from their low-lying states (at large time separation) and their 
masses will be degenerated with the masses of S- and, respectively, V^^. 

^N is the number of hypercolors. 



preserves the gauge invariance but breaks the chiral symmetry. The chiral properties 
of the theory can be restored in the continuum hmit. 

The calculation of (19) at present is not yet possible because the computation of 
detM{U) at each Monte Carlo step is very CPU time and computer memory con- 
suming. Like in QCD we assume that the quenched approximation det M{U) ~ 1 is 
a reasonable approximation also in our model. To evaluate the fermion propagator 
G{x,y) with standard Monte Carlo techniques in the quenched approximation we 
statistically integrate (19) by replacing the functional integral over the U variables 
by the measure exp{—S{U)}[dU] using the standard algorithms. 
Given a gauge field configuration U' one obtains the fermion propagator of this con- 
figuration U' by looking at the inverse fermion matrix M{U')~^ in eq. (19). The 
correlators of the different bound states will be computed by averaging eq. (16) over 
different gauge configurations. 

3.1 Matrix inversion technique 

Most of the computational effort computing the correlators goes into the construction 
of the fermion propagator G{x, y). We have to find the inverse of the fermion matrix 

M{x,y) = M{U)x,y, which obeys 

Y,M{x,y)G{y,z) = 6,,, (20) 

y 

Since M is a complex matrix of dimension]^ Vm = N^ x Nt x N x A it is not possible 
to store G{x, y) for all x and y since this involves arrays of the order of V^ complex 
numbers. Fortunately, we can exploit translation invariance such that one only needs 
to know G{x, y) for all x and for one selected point y. 
Generally, one constructs G{x,y) by solving 

J2M{x,y)G{y,z) = S,{x) (21) 

y 

where Sz{x) is some external source which can be different from 6x,z- The propagator 
G{x,y) will be the vector M~^S. 

Because we expect the fermion masses of the constituents to be much smaller than a 
typical bound energy we should invert the fermion matrix for small masses. However, 
as the fermion masses tend to the chiral limit, the matrix inversion of M becomes a 
hard problem because M tends to be singular. Therefore one usually considers a few 
heavy fermion masses and extrapolates to the chiral limit. 

In our calculation the inversion of the fermion matrix is performed by the Minimal- 
Residual algorithm or by the Conjugate Gradient algorithm in the case that the first 



^M has lattice indices, hypercolor indices and spin indices. All isospin index inversion can be 
simply done by hand to economize memory and CPU time. 



method does not converge. To economize computer memory we organize the inver- 
sion algorithms on a checkerboards even/odd-sphtting[]. 



3.2 Smearing technique 

The standard method to determine the bound state masses is to calculate the two- 
point function 

Cr{t) = Y.{0kx,t)0rm 

X 

= E{l(^|C'r(0)|0)|V<* + |(0|Oj(0)|r2)|V<(^'-*)} (22) 

n 

where Or{x) = ip{x)rilj{x) is an operator with the quantum number of the bound 
state in question, Nt denotes the number of lattice points in the time direction and 
\n) is an eigenstate of the transfer matrix. Here we consider a lattice with periodic 
boundary conditions. 

Our goal is to evaluate the mass of the lowest lying state contributing to (22). Here 
the problem of isolating this term in the sum over all eigenstates occurs. A reliable 
value for the mass can be extracted only at large time t. An efficient method to 
extract the ground state contribution to Cr is to construct an operator Or which has 
a large overlap with the ground state |no). A pointlike operator Or{x) performs very 
poorly as an operator which generates bound states. One therefore constructs a new 



bound state operator by smearing |16| the fermion field iIj{x) with a suitable wave 
function F{x,x'). One has to choose the wave function which has a physical extent 
over several lattice spacing. F has to be smooth and to vanish for large separation 






We have chosen F dependent in a gauge invariant way on the link gauge variables U. 
The smeared fermion field takes the form: 

ijsix,t)=J2F{x,x',U,t)^ix',t) (23) 

x' 

where the subscription "s" indicates smeared fields. In this work the wave function 
is chosen to be of the Gaussian type (characterized by two parameters a and n) 

F{x, x\ U, t) = (1 + aH{x, f , [/, t))" (24) 

with the hopping matrix 

H{x, x', U,t) = {U{{xx')) + U^iix'x))] |,=(i,,^),,,=(t,,^,) (25) 



^ We defined a lattice point as even (odd) if its coordinate sum {x + y + z + t) is even (odd) . Eq. 
(21) can be written in a checkerboard basis and if one knows G on one checkerboard (Ge, say) one 
can reconstruct Go- 



This wave function is defined in local time which means that it involves gauge fields 
only in one time slice and is gauge covariant. 

The smeared fermion propagator Gs{x,x') is related to the local fermion propagator 
G{x, x') by 

Gs{x, x') = J2 F\x, y, U, t) J2 G{y, z)F{z, x\ U, t) (26) 

y z 

where x = {t,x) and x' = {t',x'). We mention that the construction of the smeared 
fermion propagator does not take much more computer time than the calculation of 
the local propagator. 

The ground state dominance of the correlation functions Cr is signalled by the oc- 
currence of a plateau in the local mass 

In this case the local mass is time independent for a range of t. It was impossible 
to attain a plateau without smearing the fermion fields. In our calculation we have 
tuned the parameters a and n of the wave function (24) to obtain early plateaux in 
the local masses of the smeared operators. Defining 

^2 / j:,x'Tr[F{x,6)F^ix,0)] \ ^28) 

\ E,Tr[F{x,0)F^x,0)] / 

the optimisation of the wave function coincides in our calculation with the radius r 
approximately of 3 x a, where a is the lattice spacing. This is a resonable size for a 
bound state of a typical size less or equal the spacing a. 

3.3 Setting the scale and the hopping parameter 

From lattice calculation we can evaluate only dimensionless quantities. To connect 

these quantities to physics we have to set a scale. On the lattice we choose the scale 

to be the lattice spacing a and we express all quantities in unit of a. For small enough 

a the lattice results up to a renormalization of the parameters and the fields yield 

the continuum theory. 

In our model the experimental value of the W-boson mass will be used to setting 

the a scale. For setting the space lattice a we identify the experimental value Mw = 

80.22(26) GeV with the value of the mass of the V''^ resonance in the dimensionless 

quantity m^kj-A x a determined from the simulation. 

The lattice value of rrijkrpA x a has to be determined for heavy fermion masses to 

avoid the singularity of the fermion matrix. However, we expect that the constituent 

fermion masses are much smaller than a typical binding energy. Therefore we may 

proceed in analogy to QCD and extrapolate the quantities to the chiral limit. 

In QCD the critical value kc of the hopping parameter follows from the determination 

of the pion mass 

»' ' '"- = ^ (i - ^) '28> 

10 



suggested by the lowest order chiral perturbation theory |T^. The pion is the Gold- 
stone boson predicted by the spontaneously broken chiral symmetry and should be 
a massless particle in the chiral limit. The non- vanishing experimentally measured 
mass of the pion is the result of the breaking of flavor SU(2) symmetry. 
In our model there is no Goldstone boson which can play the role of the pion in QCD. 
However, when the chiral anomaly is switched off the pseudo scalar bound state be- 
haves like a Goldstone boson. We emphasise that swiching on the chiral anomaly by 
choosing topological non trivial gauge configurations we expect that the pseudoscalar 
bound state S"^ acquires a mass by the chiral anomaly in analogy to the r]' in QCD. 
The critical value k^ of the hopping parameter can be evaluated by assuming that 
the pseudo scalar bound state S^ behaves like a Goldstone boson and therefore it 
should be massless. In analogy to QCD one makes the following ansatz 

«' X < = ^ (^ - ^) (30) 

In our fits the linear dependence of rni^ on the inverse hopping parameter 1/k is very 
well satisfied. All other bound state masses can be evaluated by linearly extrapolating 
their lattice predictions rrar x a to the critical hopping parameter kc as suggested by 
our strong coupling expansion Jill. ^"^ particular, for the vector isotriplet V^^ mass. 



the dimensionless quantity m^kj^A x a which is used to set the scale can be extrapolated 
into the chiral regime by the linear ansatz 



1 1 

rC rCf 



a X m^krpA = B + C [- — — ) (31) 



in analogy to the linear ansatz for the p mass in QCD. 



4 The pseudoscalar isosinglet mass. 



We have shown in the ref. |]T1[ that the vector isotriplet bound state can be the 
lightest bound state provided that the pseudoscalar isosinglet acquires a mass from 
the chiral anomaly. In order that our model may be considered a viable theory of 
electroweak compositeness it has to turn out that this bound state is heavier than 
the vector isotriplet one. 

A non-perturbative estimate of the contribution of the chiral anomaly to the pseu- 
doscalar isosinglet mass is one of the most challenging problems of the lattice Monte 
Carlo simulations. In the past there were a few attempts to make a quantitative 
Monte Carlo analysis of that problem in QCD [^, |l^, ^. In this section we will 



discuss the calculation of the mass of the pseudoscalar isosinglet bound state by 
applying and improving, using the smearing technique, some methods developed in 



In QCD, in the limit of zero bare masses of the u and d quarks, the pseudoscalar 
mesons (tt and rj') are Goldstone bosons due to the spontaneous symmetry breaking 

11 



of the SU(2) axial symmetry. Because the quarks possess small but nonzero bare 
masses, the tt has a small mass and only approximates a Goldstone boson. However 
the f]' meson is too heavy for approximating a Goldstone boson. This is the U{1) 



problem |^. t'Hooft removed the puzzle by showing that the existence of topologi- 
cally non-trivial gauge configurations such as instantons resolves the f/(l) problem. 
In any confining Yang-Mills theory with fermion carrying non-trivial isospin quantum 
number the same problem arises. The pseudoscalar isosinglet bound states (the anal- 
ogous to the rj' in QCD) acquire heavy masses due to instanton effects. Their masses 
can be evaluated using the method of lattice gauge theories. There are, however, a 
number of technical difficulties which restrict the feasibility of the study. 
The propagator of the pseudoscalar isosinglet S";!^ bound state is given by equation 
(18). From the numerical evaluation of this propagator we can determine the mass 
of the S"""" bound state using the lattice spectroscopy methods presented in section 
3. The major difficulty comes from the large amount of computer time needed to 
evaluate the disconnected contribution to C°". For this purpose we need the fermion 
propagator G{x, x) for all points x , which requires N^ x Nt inversions of the fermion 
matrix using the method of the point source on a lattice of size N^ x A^^. This is a 
hopeless task because it needs to much computer time. 

The problem can be solved requiring Nf fermion matrix inversions by calculating the 
ratio 

mt) ^ ^ (32) 

for Nt walls . Here, C'^^ and C'^J^ denote the propagators computed at a suitable 
chosen value of x in eq. (17) and (18). This means that one does not perform the 
sum over x in (17) and (18). The long distance behaviour of the propagators (17-18) 
for a fixed x is then |jT^ 



c'Tit) = A[—^ + 



tD [Nt-t)^ J 

c;w - b(^^ + ^^^-^) (33) 



where A, B and D, D' are constants, m'^ and m^^^ denote the masses of S"^ and S*., 
respectively. The smearing of the source over the space-like lattice for each time wall 
has the effect to partially sum over the space-like point over each time wall and thus 
reduces the t^ and t^ dependence of the propagators (33). To a good approximation 
we expect that with smearing the exponents D and D' are Z) ^ Z^' ^ 1. The ratio 
R{t) eliminates the t^ and t^ dependence of the propagators C'^^^ and C'^^ and can 
be expressed for sufficiently large time t by 

R{t) - E (e-^™-5* + g-Am,,(JV*-t)^ (34^) 

where E = A/B is a constant and Am^^ = m^" — m^^ represents the mass difference 
between the S""" and the S^ bound states for any k < k^. Because, in the quenched 

12 



approximation, this mass difference is a pure hypergluonic effect, it does not depend 
on the mass of the Majorana fermion fields and also on the hopping parameter k. At 
the chiral limit the mass of S- is equal to zero and Am^^ correspond to the mass of 
gan^ Checking that the ratio R{t) is independent of the hopping parameter k one can 
extrapolate Am^^ to the chiral limit by a constant function in k obtaining then m'^^. 
The contribution of the disconnected term to the propagator of S*"" is present only 
if topological nontrivial gauge configurations are used. C'^J^ = C'^^^ (statistically) for 
topological trivial configurations. Topological nontrivial configurations can be ob- 
tained thermalizing an instanton on the lattice using the heat bath updating. The 
cooling algorithm is also used to determine the topological property of the configu- 
rations. Cooling p2[ works as described next: we use the standard plaquette action 
which will be locally minimized during the iterations sweeps of the heat bath algo- 
rithm. The effect is to locally smoothen the lattice gauge configuration, removing the 
ultraviolet fluctuations which are responsible for destroying the topological charge. 
The starting instanton which will be thermalized can be defined on the lattice by dis- 
cretizing the instanton solution proposed by t'Hooft [0. To control the topological 



properties of the gauge configurations we have used Peskin's lattice definition of the 



topological charge [23] in the symmetrized version [24 

Q = -E E ^^Tr[U{x),^U{x),^] (35) 

which has the right continuum limit. Here e^j^ypfj is the generalized total antisymmet- 
ric tensor^ in any direction (positive and negative) of the Euclidean space and U{x)^^ 
is the product of four links U{b) around a plaquette lying in the fiu plane and with 
starting point x. There are other definitions of the topological charge on the lattice 
20| , p5|] . We have chosen Peskin's one because it is the simplest to be programmed 



on the computer. We emphasize that a less naive (and more complicate) definition 
of the topological charge would be needed if one would like to know the topological 



charge density Q{x) or the topological susceptibility like in ||20|. However, in our 
work we only need a naive control over the global topological charge Q = Y^x£aQ{^) 
and the simple definition of Peskin is sufficient for our purpose. 

In order to perform a Monte Carlo integration with topological non-trivial configura- 
tions one has to generate a set of configurations having integer topological charges. 
Unfortunately, because the lattice definition (35) of the topological charge is a dis- 
cretization of the continuum definition which allows non integer values of Q, we can't 
generate configurations with an exact integer topological charge, but only with a Q 
lying near an integer value. We accepted only configurations with Q lying in an in- 
terval ±0.1 around an integer number. 

An other difficult is to generate configurations with |(5| > 1 because the small volume 
of the lattice does not allow a proper thermalization of such instantons: they loose 
their topological charge during the iterations and and in configurations with |Q| = 1 

®For example: 1 = £1234 = -£-1234- 
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or Q = 0. 

5 Results 

5.1 Computational details 

The simulation was performed on the Cray YMP at the ETH in Ziirich and on the 
NEC SX-3 at the CSCS in Manno. The simulation was done on different lattices in 
order to understand the finite volume and a effects (see Table 1). The SU{2) config- 
urations were generated by the combination of heat bath and over-relaxed updating 
(1 heat bath for 6 over-relaxed sweep). The cooling updating was also used for topo- 
logical non-trivial configurations to determine the topological charge. 
The plateaux of the local masses (27) were obtained from the fit to local masses of 
smeared-smeared correlators. They were independent on the parameters a and n of 
the applied smearing function (24) f\. Having established the wave function indepen- 
dence, for the rest of the calculation we performed the estimation of the spectrum 
using only the optimized wave function with a = 3 and n = 25. 
The masses were calculated at different k < k^s from a two parameter fit to eq. (22) 
or (34) on a time interval [tmin, tmax] determined by the plateaux and by studying the 
stability of the fits under changes to the fitting range. The tipical best fitting ranges 
for the different lattices A,B,C (see Table 1) are [tmin, tmax] =[6,10], [4,7] and [10,16], 
respectively. 

All needed fits were performed by minimizing the correlated x^ ^ind the statistical 
errors were calculated using the jackknife method [^ (with binning to control the 
autocorrelation) . 

5.2 Phase diagram and a^^{/3) dependence 

The critical hopping parameter kc was determined by linearly extrapolating m'i^ to 
the chiral limit as explained in section 3.3. For the different lattices the phase diagram 
in the [k, /5)-plane is plotted in Fig. 1. In this range of the parameter space there are 
no evident finite volume and a effect for these data. Data obtained from simulations 
with topological trivial and topological non-trivial configurations are consistent with 
each other. 



In a previous work [|14[ we have computed the mass m^^ by a strong coupling expan- 
sion to high orders. In Fig. 2 we plot in the {k, /3)-plane the critical line determined 
by the chiral limit m.y^{k,f3) = over a large range of (3 values (from 0.0 to 2.8). In 
the same plot we compare the critical line with the Monte Carlo data. Considering 
that the convergence radius of the strong coupling expansion can be estimated to be 
Pconv. — 1.7 the agreement between the two calculation is surprising. 



^Although the t-dependent local masses depend on the underlying wave function, they finally 
end up with the same plateau! 



14 



The lattice spacing a was determined from the vector isotriplet mass extrapolated to 
the chiral limit and normalized with the experimental value of the W-boson mass. 
In Fig. 3 we have plotted the Monte carlo data for the different lattices. We clearly 
see in the plot that the values of a~^{(3) for a fixed (3 are affected by finite volume 
effects, as expected. The prediction of the strong coupling expansion and the Monte 
Carlo data are plotted in Fig. 4. For values of (3 less than the estimated convergence 
radius Pconv. — 1-7 we obtain the function a~^(/5) (solid line). For values of /3 bigger 
than l3conv. the behaviour of this function can be interpolated by the Monte Carlo 
data (dashed line). 

5.3 Spectrum 

5.3.1 The axial vector and tensor bound state masses 

The masses of the axial vector V'^ and the tensor B^^"^ bound states are evaluated on 
the lattices Ai (i=l,...3), Bi (i=0,..3) and Ci (i=l,..3) (see Table 1) using topological 
trivial gauge configurations. After a thermalization by 5000 iterations sweeps we 
use gauge configurations separated by 200 sweeps. For 1 heat bath sweeps we do 6 
over- relaxed sweep. 

For each (3 we have fitted the masses using eq. (22) for eight different hopping 
parameters and extrapolated them to the chiral limit. 

In Fig. 5 and Fig. 6 we present the result of the lattices Bl, B2, B3 and CI, C2, C3 
for which the finite volume effects are small. The same results are presented also in 
Table 2. We plot the masses as a function of the lattice spacing a normalized with 
the W-boson mass at 80 GeV. All data points with itiy x a > 1 are excluded from the 
analysis of the finite a effects. These correspond to masses evaluated on the lattice 
BO with the smallest j3. The continuum limit of the masses is evaluated by linearly 
extrapolating to a = for the different lattices. The results of the extrapolations 
are presented in Fig. 5 and Fig. 6. Finite volume effects at the continuum limit 
are negligible between the lattices Bl, B2, B3 and CI, C2, C3 in comparison to the 
corresponding statistical errors. 

5.3.2 The pseudoscalar isosinglet mass 

The mass of the pseudoscalar isosinglet S*"" is evaluated on the lattices BO*, Bl* and 
B2*. We used topological non-trivial gauge configurations with topological charge 
Q = 1. After a thermalization by 5000 iterations sweeps we used gauge configura- 
tions separated by 1000 sweeps. The evolution of the topological charge during the 
iterations is plotted in Fig. 7. 

For each (3 we have fitted the mass difference Atti^^ with the ratio R{t) of eq. (34) for 
five different hopping parameters k. In Fig. 8 we plot the different values of Atti^^ x a 
as a function of the m^^ x a^ (which is proportional to the inverse hopping parame- 
ter). The same result is presented also in Table 3. We note that it is independent of 
mi^ , as it should be. After checking this mass independence the extrapolation to the 
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chiral limit at rn^ = was done by a constant function. 

In Fig. 9 we plot the chiral limit of Am^g as a function of the lattice spacing a. 
At the chiral limit the mass difference coincides with the mass of the pseudoscalar 
isosinglet m"". We clearly see that decreasing the lattice spacing a the pseudoscalar 
isosinglet mass increase and it is bigger than the W-boson mass. 

6 Conclusion 

We have discussed a composite model for the weak bosons of the SM. The model is 
based on a Yang-Mills theory without Goldstone bosons. In this theory the W-boson 
of the SM is represented by the vector isotriplet bound state. As a main result our 
lattice Monte Carlo simulation has shown that the vector isotriplet bound state is 
the lightest bound state in our model as it should be for any viable candidate of 
electroweak composite model. 

We have developed an efficient method to compute the mass of the pseudoscalar 
isosinglet bound states. This computation requires the evaluation of disconnected 
fermion loops and the generation of topological non-trivial gauge configurations. 
We have also predicted the mass of the first bound states heavier than the vector 
isotriplet. These bound states are a vector isotriplet and a vector isosinglet and a 
pseudoscalar with masses in the range of a few hundred GeV (see Fig. 5, Fig. 6 and 
Fig. 9). These predictions open new experimental perspective at LEPII and LHC 
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Figure Caption 

1. Lattice Monte Carlo data for the phase diagram in the {f3, A;)-plane. 

2. Lattice Monte Carlo data and strong coupling expansion prediction for the 
phase diagram in the {(3, /c)-plane. The dotted line represents the limit of 1/8 

for /9 — s> oo. 

3. Lattice Monte Carlo data for the inverse lattice spacing a~^ as a function of f3. 

4. Lattice Monte Carlo data and strong coupling expansion calculation for the 
inverse lattice space a^^ as a function of (3. The dotted line represents the 
estimate convergence radius of the strong coupling expansion. The dashed line 
is a qualitative estimation of a~^{(3) based on an interpolation of the Monte 
Carlo data of the biggest lattice. 

5. Lattice Monte Carlo prediction of the triplet BA. bound state mass as a function 
of the lattice spacing. Dashed and dotted lines are the extrapolation to the 
continuum limit. 

6. Lattice Monte Carlo prediction of the axial vector V^ bound state mass as a 
function of the lattice spacing. Dashed and dotted lines are the extrapolation 
to the continuum limit. 

7. Evolution of the topological charge during the iterations. The topological charge 
is measured at each 50 iterations. 

8. The mass difference Am^" as a function of the m'i mass in unit of the lattice 
spacing a. The dashed lines represent the extrapolations to the chiral limit. 

9. The pseudoscalar isosinglet mass as a function of the lattice spacing a. 
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lattice 


/9 


Ns, Nt 


no. config's. 


a-^GeV 


kc 


Al 


2.3 


6 


24 


44 


120(10) 


0.164(6) 


A2 


2.4 


6 


24 


40 


145(12) 


0.156(6) 


A3 


2.5 


6 


24 


40 


156(24) 


0.153(4) 


BO* 


2.2 


8 


16 


40 


129(7) 


0.173(1) 


BO 


2.2 


8 


16 


40 


119(9) 


0.173(1) 


Bl* 


2.3 


8 


16 


40 


156(4) 


0.170(2) 


Bl 


2.3 


8 


16 


40 


141(8) 


0.166(2) 


B2* 


2.4 


8 


16 


40 


181(17) 


0.154(2) 


B2 


2.4 


8 


16 


42 


174(16) 


0.154(1) 


B3 


2.5 


8 


16 


40 


198(17) 


0.153(2) 


B4 


2.7 


8 


16 


40 


******* 


0.152(8) 


CI 


2.3 


1 


2, 36 


40 


158(6) 


0.1672(4) 


C2 


2.4 


1 


2, 36 


40 


177(4) 


0.1558(3) 


C3 


2.5 


1 


2, 36 


40 


225(12) 


0.1515(4) 



Table 1: Parameters of the lattices used for this work. The inverse lattice spacing 
is obtained by fixing the vector isotriplet mass to 80 GeV. The critical hopping 
parameter is obtained by extrapolating m^^ to the chiral limit. Asterisks (*****) 
indicate that the fit was not accepted due to a large x^ ot the missing of a plateau 
in the local masses in the region of the fit. The lattices BO*, Bl* and B2* indicate 
simulations performed with topological non-trivial configurations. 



lattice 


a-VGeV 


m^jkrpi GeV 


m^j.^fe GeV 


Al 
A2 
A3 


120(10) 
145(12) 
156(24) 


73(47) 
76(21) 
111(40) 


83(48) 

104(27) 

98(39) 


BO 
Bl 
B2 
B3 


119(9) 
141(8) 
174(16) 
198(17) 


202(36) 
108(27) 
135(25) 
122(26) 


152(37) 
125(24) 
130(31) 
128(29) 


CI 
C2 
C3 


158(6) 

177(4) 

225(12) 


97(14) 
96(10) 
113(20) 


128(27) 
163(34) 
133(20) 



Table 2: Lattice Monte Carlo predictions of the triplet BL and singlet V^ bound 
states as a function of the inverse lattice spacing a~^. 
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lattice 


a-VGeV 


ml^^ X a^ 


Am";" X a 


BO* 


129(7) 


0.270(27) 


0.71(15) 






0.359(20) 


0.71(14) 






0.442(18) 


0.70(14) 






0.562(16) 


0.69(12) 






0.719(17) 


0.68(11) 


Bl* 


156(4) 


0.489(6) 


0.98(15) 






0.543(7) 


1.00(16) 






0.596(7) 


1.02(17) 






0.632(8) 


1.03(17) 






0.702(9) 


1.05(18) 


B2* 


181(17) 


0.535(14) 


0.97(9) 






0.610(14) 


1.00(10) 






0.712(13) 


1.06(12) 






0.776(12) 


1.10(13) 






0.830(11) 


1.14(14) 



Table 3: Lattice Monte Carlo predictions of the mass difference Am"" in unit of the 
lattice spacing. 
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